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Abstract 

^ ' We calculate the cross sections of inclusive J/ ip production in photoproduction 

and two-photon scattering, involving both direct and resolved photons, and in e~^e 
Qh' annihilation at next-to-leading order (NLO) within the factorization formalism of 

nonrelativistic quantum chromodynamics (NRQCD), including the full relativistic 
corrections due to the intermediate ^S^\ ^Sf\ and 3pjsi color-octet (CO) states. 
Exploiting also our previous results on hadroproduction, we perform a combined fit 
J> . of the respective CO long-distance matrix elements (LDMEs) to all available high- 

quality data of inclusive J/^p production, from KEKB, LEP II, RHIC, HERA, the 
QQ ' Tevatron, and the LHC, comprising a total of 194 data points from 26 data sets. 

^ ' PACS numbers: 12.38.Bx, 13.60.Le, 13.85.Ni, 14.40. Gx 
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The factorization formalism of NRQCD provides a rigorous theoretical framework 
for the description of heavy- quarkonium production and decay. This implies a separa- 
tion of process-dependent short-distance coefficients, to be calculated perturbatively as 
expansions in the strong-coupling constant a^, from supposedly universal LDMEs, to be 
extracted from experiment. The relative importance of the latter can be estimated by 
means of velocity scaling rules; i.e., the LDMEs are predicted to scale with a definite 
power of the heavy-quark {Q) velocity v in the limit f ^ 1. In this way, the theoretical 
predictions are organized as double expansions in and v. A crucial feature of this for- 
malism is that it takes into account the complete structure of the QQ Fock space, which 
is spanned by the states n = 25+i2^W ^[\;[^ definite spin S, orbital angular momentum L, 
total angular momentum J, and color multiplicity a = 1,8. In particular, this formalism 
predicts the existence of CO processes in nature. This means that QQ pairs are produced 
at short distances in CO states and subsequently evolve into physical, color-singlet (CS) 
quarkonia by the nonperturbative emission of soft gluons. In the limit f — ?• 0, the tra- 
ditional CS model (CSM) is recovered in the case of S'-wave quarkonia. In the case of 
J/iIj production, the CSM prediction is based just on the ^5*^ CS state, while the leading 
relativistic corrections, of relative order 0{v^), are built up by the ^S^o \ ^S? , and ^P^ 
(J = 0, 1, 2) CO states. The CSM is not a complete theory, as may be understood by 
noticing that the NLO treatment of P-wave quarkonia is plagued by uncanceled infrared 
singularities, which are, however, properly removed in NRQCD. 

The test of NRQCD factorization has been identified to be among the most exigent 
milestones on the roadmap of quarkonium physics at the present time [2]. While, for 
J/ip polarization, comparisons of HERA and Tevatron data with NRQCD predictions, 
which are not yet fully known at NLO, unravel a rather confusing pattern, the situation 
is eventually clearing up for the J/ip yield, which is now fully known at NLO in NRQCD 
for direct photoproduction [3] and hadroproduction [IfS] . In fact, it has been demon- 
strated [1] that the set of CO LDMEs fitted to transverse-momentum (pt) distributions 
measured at HERA [GfTj and by CDF at Tevatron II [H] also lead to very good descrip- 
tions of distributions in the jp cm. energy W and the inelasticity z, which measures the 
fraction of 7 energy passed on to the J/ip meson in the p rest frame, from HERA piT] 
and of pt distributions from RHIC [9] and the LHC [10] . On the other hand, the Teva- 
tron II [S] data alone can only pin down two linear combinations of the three CO LDMEs 
[SfTT] . and the fit results of Ref. |5j are incompatible with Ref. It is the purpose 
of this Letter, to overcome this highly unsatisfactory situation jeopardizing the success 
of NRQCD factorization by performing a global fit to all available high-quality data of 
inclusive unpolarized J/ ip production, comprising a total of 194 data points from 26 data 
sets. Specifically, these include pt distributions in hadroproduction from PHENIX |9j 
at RHIC, CDF at Tevatron I [12] and Tevatron II [g], ATLAS [13j, CMS [ID], ALICE 
[13], and LHCb [TB] at the LHC; p"^, W, and z distributions in photoproduction from 
ZEUS [16J and HI [6] at HERA I and HI [7] at HERA II; a p'^ distribution in two-photon 
scattering from DELPHI [17j at LEP II; and a total cross section in e~^e~ annihilation 
from Belle [18] at KEKB. 

Incoming photons participate in the hard scattering either directly or via partons into 
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{O-^/^i^S^o)) I (4.97 ± 0.44) X 10-2 GeV^ 
(2.24 ± 0.59) X 10-=^ GeV^ 

(C)W(3pM)) (_1,51 ± 0.20) X 10^2 Q^y5 



Table 1: NLO fit results for the J ftp CO LDMEs. 



which they fiuctuate (resolve) intermittently, and both modes of interaction contribute 
at the same order of perturbation theory. Therefore, we need to extend the theoretical 
ingredients available from Refs. [3lH] by also treating 7P — )■ J/ip + X with the photon 
being resolved and 77 — )■ J/ip + X with none [19], one, or both of the photons being 
resolved at NLO in NRQCD. We repeat the analysis of Ref. [19], in which the Coulomb 
singularities were regularized by v, using dimensional regularization as in Refs. [3]!^ in 
order to obtain analytic expressions sufficiently compact for our purposes. We also find it 
necessary to revisit e~^e~ J/ip + X ai NLO in NRQCD because the results of Ref. ^20j 
have not yet been verified by an independent calculation, are only available in numerical 
form, and lack the ^S'f contribution, which comes both with X = qq |2l] and gg. Higher- 
order corrections to the CSM process e~^e~ — )■ cc[^S'f 'j^f^f |22], which enters our analysis 
at NLO, are beyond the order considered here. 

The additional analytic calculations proceed along the lines of Refs. [314] and are not 
described here in detail for lack of space. We merely present our master formula based 
on the factorization theorems of the QCD parton model and NRQCD [1]: 

da{AB J/ip + X) = j dxidx2dyidy2 fi/A{.Xi)fk/i{yi)fj/B{.X2)fi/j{y2) 

i,j,k,l,n 

X (C'^/^[n])rfa(A;/ ^ cc[n] + X), (1) 

where fi/A{xi) is the parton distribution function (PDF) of parton i = g,q,q in hadron 
A = p,p or the fiux function of photon i = 7 in charged lepton A = e~,e^, fk/iivi) is 
SikS{l — yi) or the PDF of parton k in the resolved photon i, da{kl — )■ cc[n] + X) are the 
partonic cross sections, and {0'^^'^[n\) are the LDMEs. In the fixed-fiavor- number scheme, 
we have q = u,d, s. In the case of e^e~ annihilation, all distribution functions in Eq. 
are delta functions. As in Refs. [3]H] . X always contains one hard parton at leading order 
(LO) and is void of heavy fiavors, which may be tagged and vetoed experimentally. 

We now describe our theoretical input for our numerical analyses. We set rric = 
1.5 GeV, adopt the values of rrig, a, and the branching ratios B{J/iIj e^e~) and 
B{J/ip — 7- yU^/i^) from Ref. |23], and use the one-loop (two-loop) formula for a^J^^\fi), with 
Uf = A active quark fiavors, at LO (NLO). As for the proton PDFs, we use set CTEQ6L1 
(CTEQ6M) [23] at LO (NLO), which comes with an asymptotic scale parameter of Aq^j-, = 
215 MeV (326 MeV). As for the photon PDFs, we employ the best-fit set AFG04_BF of 
Ref. [25]. We evaluate the photon fiux function using Eq. (5) of Ref. [26], with the 
upper cutoff on the photon virtuality chosen as in the considered data set. As for the 
CS LDME, we adopt the value (O-'/^i^sP)) = 1.32 GeV^ from Ref. [27]. Our default 
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choices for the renormahzation, factorization, and NRQCD scales are Hr = Hf = rriT and 

/iA = rUc, respectively, where mx = \J'Pt ^ is the J/tp transverse mass. The bulk of 
the theoretical uncertainty is due to the lack of knowledge of corrections beyond NLO, 
which are estimated by varying /x^, /ij, and /ia by a factor 2 up and down relative to their 
default values. 

We exclude from our fit all data points of photoproduction and two-photon scattering 
with pt <1 GeV and of hadroproduction with < 3 GeV, which cannot be successfully 
described by our fixed-order calculations as expected. This leaves a total of 194 data 
points. The fit results for the CO LDMEs obtained at NLO in NRQCD with default scale 
choices are collected in Table [H They depend only feebly on the precise locations of the 
Pt cuts. In the following, we use the values of Table [H throughout. 

In Figs. [Hand [2]^a), all data sets fitted to, except the single data point from Belle [18], 
are compared with our default NLO NRQCD results (solid lines). For comparison, also the 
default predictions at LO (dashed lines) as well as those of the CSM at NLO (dot-dashed 
lines) and LO (dotted lines) are shown. The yellow and blue (shaded) bands indicate the 
theoretical errors on the NLO NRQCD and CSM predictions. We observe from Figs. [T] 
that the experimental data are nicely described by NLO NRQCD, being almost exclusively 
contained within its error bands, while they overshoot the NLO CSM predictions typically 
by 1-2 orders of magnitude for hadroproduction and a factor of 3-5 for photoproduction. 
The description of the z distributions in photoproduction by NLO NRQCD significantly 
benefits from two features, rendering it considerably more favorable than in Refs. [3)i4] . On 
the one hand, as conjectured in Refs. [3P] . resolved photoproduction usefully enhances the 
cross section in the low- 2; range, being dominant for 2^0.25, as is evident from Fig. [2](b). 
On the other hand, owing to the negative value of {O-'/^ {^Pf^)) in Tabled the ^Sf and 
^Pf^ contributions interfere destructively thus attenuating the familiar rise in cross section 
in the limit 2; — ?• 1, as may be seen from Fig. [2](c). As for the p\ and W distributions in 
photoproduction, the cut z > 0.3 (0.4) applied by HI (ZEUS) greatly suppresses resolved 
photoproduction, to the level of 1%. In contrast to the LO analysis of Ref. [2B], the 
DELPHI [17] data tend to systematically overshoot the NLO NRQCD result, albeit the 
deviation is by no means significant in view of the sizeable experimental errors. As is 
evident from Fig. [2](d), this may be attributed to the destructive interference of the '^S^'^ 
and ^Pf' contributions mentioned above, which is a genuine NLO phenomenon. We have 
to bear in mind, however, that the DELPHI measurement comprises only 16 events with 
Pt > I GeV and has not been confirmed by any of the other three LEP II experiments. In 
two-photon scattering at LEP II, the single-resolved contribution vastly dominates over 
the direct and double-resolved ones, as was already observed for the LO case in Ref. [28] . 
The Belle measurement, (T(e+e~ — )■ J/ip + X) = (0.43 ±0.13) pb, is compatible both with 
the NLO NRQCD and CSM results, (O.TO+J;^^) pb and (0.24+[|;J°) pb, respectively; at 
LO, where X = g, we are dealing with a pure CO process, with a total cross section of 
0.23 pb. The overall goodness xlo.i. = 857/194 = 4.42 of our NLO NRQCD fit, which we 
quote for completeness, is of limited informative value, since the theoretical uncertainties 
exceed most of the experimental errors. 
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Figure 1: NLO NRQCD fit compared to RHIC ^, Tevatron [gllT^ . LHC 
and HERA piTllTB] data. 
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Figure 2: (a) NLO NRQCD fit compared to LEP II ^17j data and (d) its decomposition 
to cc[n] channels. Decomposition of z distribution at HERA I |I6] (b) to contributions 
due to direct and resolved photoproduction and (c) to cc[n\ channels. 



Our theoretical predictions refer to direct J /ip production, as the data from Tevatron I 
[12] do, while the data from KEKB [18], Tevatron II [8], and LHC pOTlMillTS] comprise 
prompt events and those from LEP II [T7] , HERA [6|7|IT6] , and RHIC [9] even non-prompt 
ones. However, the resulting error is small against our theoretical uncertainties and has 
no effect on our conclusions. In fact, the fraction of J/ip events originating from the 
feed-down of heavier charmonia only amounts to about 36% [T2| for hadroproduction, 
15% [7j for photoproduction at HERA, 9% for two-photon scattering at LEP II [19], 
and 26% for e^e~ annihilation at KEKB [22], and the fraction of J/ip events from B 
decays is neghgible RHIC, HERA [7], and LEP II [19] energies. Refitting the data with 
the estimated feed-down contribution subtracted yields {O'^^^i^S)^^)) = (3.04 ± 0.35) x 
10-2 GeV^ (O-^/^i^sP)) = (1.68 ± 0.46) x 10-^ GeV^, and (^^/^(^pj®])) = (-9.08 ± 
1.61) X 10-3 GeV^ with a slightly reduced value xio.f. = 725/194 = 3.74. 

In conclusion, we performed a NLO NRQCD analysis of all available high-quality data 
of inclusive unpolarized J/ip production, from KEKB pLSj, LEP II [ITj, RHIC [9j, HERA I 
[6IT6] and II [7], Tevatron I [12] and II [8], and the LHC [T0|13|14p5] . comprising a total 
of 194 data points from 26 data sets. The fit values of the CO LDMEs in Table [1] agree 
with our previous ones [1], extracted just from the pt distributions of Refs. [6|7|8] . within 
the errors of the latter, but the new errors are about 40% smaller. In compliance with the 
velocity scaling rules of NRQCD [1], these values are approximately of order 0{v'^) relative 
to (^'^/'^(^S'f' )). This manifestly consolidates the verification of NRQCD factorizatization 
for charmonium and provides rigorous evidence for LDME universality and the existence 
of CO processes in nature. 

This work was supported in part by BMBF Grant No. 05H09GUE, DFG Grant No. 
KN 365/6-1, and HGF Grant No. HA 101. 
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